function Fig2_b()
lb = 25.66 / 6^0.5;
EC = 56.1 * 6^0.5;

% This proccessed data is processed by defect_A.ipynb
% defect_A.ipynb is located at .../processed data
sim = readmatrix('processed data/result_simZLL.csv');
sim_r = linspace(0,6,201)' * lb;
sim_e = -40:0.1:10;
sim_data = sim';

% This proccessed data is processed by defect_A.ipynb
% defect_A.ipynb is located at .../processed data
data = readmatrix('processed data/defect_A_map.csv');
exp_r = -1*data(2:end,1);
exp_e = data(1,2:end) * 1000 +110;
exp_data = ( rot90(data(2:end,2:end),3) );

% This proccessed data is processed by defect_A.ipynb
% defect_A.ipynb is located at .../processed data
curve = readmatrix('processed data/result_simPotential.csv');
curve = curve * EC;

plot_limit = [-50,50,-35,0];

data = 48;
Fcolor='k';
Bcolor='w';
Fsize=18;  %font size

figure()
ax1 = axes;
imagesc(ax1,sim_r,sim_e,sim_data)
axis(ax1,plot_limit)

xl=xlabel('position (nm)');yl=ylabel('energy (meV)');
ax.FontSize=Fsize;
yl.FontSize=Fsize;
xl.FontSize=Fsize;

ax2 =axes;
imagesc(ax2,exp_r,exp_e,exp_data)
axis(ax2,plot_limit)

ax1.YDir='normal';ax1.XDir='normal';
ax2.YDir='normal';ax2.XDir='normal';


ax.CLim=[0,3];
cMap=getColor('RdBu');
colormap(ax2,cMap);

% ax1.Visible = 'Off';
% ax1.XTick = [];
% ax1.YTick = [];
set(ax2,'color','none','visible','off');


% figure()
hold on
plot([-1*flipud(sim_r);sim_r],[flipud(curve);curve],'LineWidth',2)
axis(ax2,plot_limit)


end